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Abstract 

We consider the design problem for a Marx generator electrical network, a pulsed 
power generator. The engineering specification of the design is that a suitable res- 
onance condition is satisfied by the circuit so that the energy initially stored in a 
number of stage capacitors is transferred in finite time to a single load capacitor 
which can then store the total energy and deliver the pulse. We show that the com- 
ponents design can be conveniently cast as a structured real eigenvalue assignment 
with significantly lower dimension than the state size of the Marx circuit. Then we 
comment on the nontrivial nature of this structured real eigenvalue assignment prob- 
lem and present two possible approaches to determine its solutions. A first symbolic 
approach consists of the use of Grobner basis representations and allows to compute 
all the (finitely many) solutions. A second approach is based on convexification of 
a nonconvex optimization problem with polynomial constraints. We show that the 
symbolic method easily provides solutions for networks up to six stages while the 
numerical method can reach up to seven and eight stages. We also comment on the 
conjecture that for any number of stages the problem has finitely many solutions, 
which is a necessary assumption for the proposed methods to converge. We regard 
the proof of this conjecture as an interesting challenge of general interest in the real 
algebraic geometry field. 



1 Introduction 

Electrical pulsed power generators have been studied from the 1920s with the goal to 
provide high power electrical pulses by way of suitable electrical schemes that are slowly 
charged and then, typically by the action of switches, are rapidly discharged to provide 
a high voltage impulse or a spark (see, e.g., [S]). Marx generators (see O §3.2] or [7] 
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for an overview) were originally described by E. Marx in 1924 and correspond to circuits 
enabling generation of high voltage from lower voltage sources. While many schemes 
have been proposed over the years for Marx generators, a recent understanding of certain 
compact Marx generators structures [6] reveals that their essential behavior can be well 
described by a suitable LC ladder network where certain components should be designed 
in order to guarantee a suitable resonance condition. In turn, such a resonance condition 
is known to lead to a desirable energy transfer throughout the circuit and effective voltage 
multiplication which can then be used for pulsed power generation. 

In this paper we follow-up on the preliminary results of [H [3S], which are strongly mo- 
tivated by the observations in [S], and we transform a compact Marx generator design 
problem into a structured real eigenvalue assignment problem for a linear system associ- 
ated with the Marx network and only depending on the number of its stages. The problem 
of real pole (or eigenvalue) assignment for linear systems has been largely studied in the 
1990s, see [SB] for a survey. In its simplest form, it can be stated as follows: given ma- 
trices A e M"^", B e C e W^"^ and a monic polynomial q{s) G R[s] of degree 
n, find an m x p real matrix F such that det (sJ„ — A — BFC) = q{s) where J„ denotes 
the identity matrix of size n. This problem has generically a solution if mp > n, and it 
has generically no solution if mp < n. The situation mp = n is much more subtle. For 
this situation, it was proved in [T^ [TB] that the pole placement map from the feedback 
matrix F to the characteristic polynomial q{s) is generically not surjective. It means that 
there is a non-empty open subset of real matrices A,B,C for which there exist open sets 
of pole configurations symmetric w.r.t. the real axis which cannot be assigned by any real 
feedback. In this paper, we do not consider static output feedback pole assignment in the 
form just described, but in a structured form. The number of degrees of freedom (number 
of entries in the feedback matrix) is equal to n, the number of poles to be assigned, so it 
bears similarity with the difficult case mp = n of static output feedback pole assignment 
described above. 

Within the context of the above cited works \T2\ [T5] , the structured pole assignment 
problem that we characterize in this paper corresponds to the case where C = /„ and F is 
diagonal, which is a case not covered by these works. For this problem, we provide in this 
paper a few general results characterizing its solutions, and then show how one can tackle 
the solution. For computationally tractable first technique, based on Grobner 

basis representations, allows determining the full set of solutions, or in other words all 
the possible selections of the circuit components leading to the desired energy transfer. 
A second numerical technique, based on nonconvex polynomial optimization tools, allows 
determining one solution which is optimal in the sense that it minimizes a certain poly- 
nomial cost function. In the paper we state assumptions ensuring that the solution set is 
nonempty (that is, there exists at least one choice of the parameters leading to the desired 
resonance condition) and that the number of solutions is actually finite, namely the set of 
all solutions is a zero-dimensional set. These assumptions imply desirable termination and 
convergence properties of the symbolic and numerical techniques which are highlighted in 
the paper. Interestingly, numerical evidence reveals that these assumptions are satisfied 
for all the numerically tractable cases that have been considered, however their general 
validity is an interesting open problem. 

Additional work related to the techniques reported here corresponds to [S] and references 
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therein, where algebraic techniques (ehmination theory using resultants) are used in the 
design of power electronic devices. In this reference the authors show that parasitic 
higher order harmonics in a multilevel converter can be removed by solving a structured 
polynomial system of equations featuring a high degree of symmetry. 

The paper is structured as follows: in Section [2] we described the Marx generator network, 
we formalize the design problem and we provide the main theorem of the paper, trans- 
forming the design into an eigenvalue assignment problem. In Section [3] we present the 
symbolic solution to the eigenvalue assignment problem based on Grobner basis represen- 
tations. In Section m we present the numerical solution based on convexified formulations. 
The proof of the main theorem is reported in Section \5\ at the end of the paper for a sim- 
plified reading. Finally, Section |6] contains some discussions about the above mentioned 
conjectures and the relevance of this study as a challenge within the algebraic geometric 
field. 

Notation: Given a square matrix A, cr{A) denotes its spectrum, namely the set of its 
complex eigenvalues. Given a vector / G M", Q[/] denotes the set of all polynomials with 
rational coefficients in the indeterminates /. 



2 Marx generator design 

2.1 Circuit description and problem statement 

We consider the Marx generator network shown in Figure [1] consisting in n stages (and 
77, + 1 loops) where, disregarding the two rightmost components of the figure, each one 
of the n stages consists in 1) an upper branch with a capacitor and an inductor and 
2) a vertical branch with a capacitor only. Following [H [33 E], we assume that all the 
capacitors and inductors appearing in the upper branches are the same (corresponding to 
some fixed positive reals c and i). We will call these capacitors "storage capacitors" in the 
sequel, for reasons that will become clear next. The design problem addressed here is the 
selection of the vertical capacitors, which are exactly n, where n is the number of stages of 
the Marx circuit. We will call these capacitors "parasitic capacitors" due to their position 
resembling that of parasitic capacitors in transmission lines. Despite their name, the 
parasitic capacitors q, i = 1, . . . ,n are not necessarily arising from any parasitic effects 
and their values will be selected in such a way to ensure a suitable resonance condition 
as clarified next. 



q 



t^c(n-l) 
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q 
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zc/n 



Figure 1: The passive circuit used as a Marx generator. 
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Following [111351 16], the inductance and capacitor appearing in the rightmost loop take the 
values n(. and c/n, respectively. We call this capacitor the "load capacitor". This selection 
preserves the resonance property (so that the product of any adjacent capacitor /inductor 
pairs is always ic) in addition to ensuring that the load capacitor is n times larger than 
each one of the storage capacitors. The problem addressed in this paper (resembling that 
one tackled in [H [351 E|) is the following. 

Problem 1 Consider the circuit in FigureUlfor given n and certain values of c and i. 
Select positive values Ci > 0, i = 1, . . .n of the parasitic capacitors and a time T > 
such that, initializing at t = all the storage capacitors with the same voltage v{0) = Vo 
and starting from zero current in all the inductors and zero voltage across the parasitic 
capacitors and the load capacitor, the circuit response is such that at t = T all voltages 
and currents are zero except for the voltage across the load capacitor. 

Remark 1 The resonance condition required in Problem [1] is a key feature to allow the 
use of the circuit of Figure [T] for pulsed power generation. In particular, assuming that 
a number of switches are used to open the loops of the circuit, the n storage capacitors 
can be charged at the same voltage Vo- Due to the property in Problem [H after T times, 
all the energy initially stored in the storage capacitors will be concentrated in the load 
capacitor. In particular, since this is a lossless circuit, we will have 

1=1 

which clearly implies vl{T) = nvo, namely both the voltage and energy transferred to the 
load capacitor is n times larger than the voltage and energy initially stored in each one 
of the storage capacitors. o 



2.2 Solution via structured eigenvalue assignment 

In this section we show that a solution to Problem [1] can be determined from the solution 
of a suitable structured eigenvalue assignment problem involving a matrix B G M"^" 
defined as 



B :-- 
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and an arbitrary set of even harmonics of the fundamental frequency luq = ^/JJc}^ to 
be assigned to the circuit. The following is the main result of this paper, whose proof is 
given in Section [51 
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Theorem 1 Consider any set of n distinct positive even integers a = (ai, . . . the 
matrix B in (QP and any positive definite real diagonal solution F = diag(/i, . . . , /„) to 
the structured eigenvalue assignment problem 

a{BF) = {al-l,...,al-l}. (2) 

Then for any value of the selection Ci = c/ fi, i = l,...^n, solves Problem U\ with 
^ ~ VTc- 

Theorem [1] shows that a solution to Problem [1] can be determined by solving an eigen- 
value assignment problem with decentralized feedback (because matrix F is diagonal). A 
generalized version of this structured pole assignment problem was studied in [31] (in- 
deed, using the notations there, the pole assignment problem in ([2]) is obtained by setting 
r = n, C = /„, mj = Pi = 1, i = 1, . . . , n). It is shown in [3l] that generic pole assignment 
depends on the dimension of a product Grassmannian, see [Ml Equ. (17)]. In our case 
it is equal to n! which is always even, and from [SU Theorem 4.2] it follows that generic 
pole assignment cannot be achieved. A geometric condition that ensures generic pole 
assignment is given in [3H Prop. 4.2] but we do not know whether this condition can be 
checked computationally. In any case it is an evidence that the question of existence of a 
real solution to our inverse eigenvalue problem does not appear to be trivial. 

Remark 2 Bearing in mind that F is diagonal, real and positive definite, the inverse 
eigenvalue problem can be equivalently cast as a symmetric inverse Jacobi eigenvalue 
problem (IJEP) (see, e.g., [HI [10]) by performing the coordinate change with T = a/F, 
which leads to 

M := T{BF)T'^ = \^BVF, (3) 

where it is evident that M = > 0. In particular, positive definiteness of M arises 
from positive definiteness of B which is established in the proof of Lemma [31 Since Q is a 
change of coordinates, then assigning the spectrum {af — 1, . . . , a,^ — 1} to the matrix M is 
equivalent to solving the eigenvalue assignment problem with the possible advantage 
that M is a symmetric matrix. o 

The following result shows that any solution to is physically implementable as it 
corresponds to positive values of the parasitic capacitors. 

Lemma 1 Any solution to the structured eigenvalue assignment problem ^ in Theo- 
remUl is such that F > 0. 

Proof: Given any solution to i^, all the diagonal entries of F are necessarily nonzero 
otherwise BF would be rank deficient which contradicts the fact that a{BF) only has 
strictly positive elements. Define F = FD, where F > is diagonal and D is a diagonal 
matrix whose elements are either 1 or —1. 

Assume now that the statement of the proposition is not true so that D has at least one 
negative entry. We consider the coordinate transformation \/F BFVf = Vf B\fF D =: 
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MD, where the last inequahty follows from the fact that F and D are both diagonal. Due 
to the coordinate transformation, all the eigenvalues of MD are positive (because they are 
the same as those of BF). Since -B > (see the proof of Lemma [3] in Section we have 
that M := VfB^/F is symmetric positive definite, so that there exists an orthogonal 
matrix Q such that M = QAQ^ , with A diagonal positive definite. Then it follows that 
MD = AQ^D which, pre- multiplied by (-\/A)^^ and post-multiplied by Q^/A leads to 

s := {Va)-^q^mdq\^ = Vaq'^dqVa, 

which establishes that the matrix S is symmetric (rightmost term) and has positive eigen- 
values (middle term which is a coordinate transformation from BF). Nevertheless, the 
matrix D is not positive definite by assumption, which leads to a contradiction. □ 



2.3 Two equivalent eigenvalue assignment formulations 



Selecting the diagonal entries / = [fi - ■ ■ fn]^ to solve the eigenvalue assignment ([2]) 
amounts to solving a finite set of n equations (polynomial in the unknown / with rational 
coefficients), each of them corresponding to one coefficient of the following polynomial 
identity in the variable s: 



n 

det{sI-BF) = n (' 



[a, 



Vs G 



(4) 



As an example, for the case n = 2, where BF 



2/1 -fl 

-Si I/2 



if one makes the simple 



selection a = (2,4), then the set of polynomial equations ensuring (jl]) corresponds to: 



2/1/2 

2/1 + t/2 



(a? 



I)(a2 - 1) = 45 
1) + (a^o - 1) = 1^ 



(5) 



In the general case, for a fixed value of n and fixed values in a, one can write a system of 
n polynomial equations in the variable / with rational coefficients, namely 



P.(/) = 0, ^ 



n. 



(6) 



In this paper, we will also adopt an alternative formulation of the problem which appears 
to be more suitable for the numerical optimization techniques developed in Section m 
The alternative formulation corresponds to inverting the eigenvalue assignment problem 
([2]), thereby obtaining an alternative set of polynomial equations in the unknowns k = 
[ki - ■ ■ /c„,]^ = [/f ^ ■ ■ ■ fn^f" with rational coefficients, which have the advantage of being 



linear in the capacitor values, indeed, ki = Ci/c, i 
problem, equation (jl]) becomes 



n. In particular, for the inverse 



n 

det(s/ - KB~^) = n - ("*^ ~ 1)^0 ' ^ 

i=l 



(7) 



where K = diag(A;) = F ^ is a diagonal matrix whose diagonal elements are the scalars 



to be determined. Similar to above, for the case n = 2, we have KB 



k2 



which, 
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for the simple selection a = (2, 4) leads to the following set of polynomial inequalities 
ensuring ([7]) (and, equivalently, (jl])): 



fa? 



1)' 
1)' 



j_ 

45 



In the general case, for a fixed value of n and a fixed values in a, one can write a system 
of n polynomial equations in the variable k with rational coefficients, namely 



qi{k) = 0, i 



n. 



(9) 



Formulation ([7]), will be used in Section H] due to the advantageous property that all 
entries of each solution to this polynomial system are in the interval (0, 1) as established 
in the next lemma. 

Lemma 2 Given anyn > 1 and any set of distinct positive even integers a = (ai, . . . , 
each solution k = {ki, . . . , kn) to the inverse eigenvalue problem ^ is such that < fcj < 1 
for all i = 1, . . . ,n. 



Proof: The fact that ki > for all i = l,...,n follows from Lemma [T] and K = 
diag{ki,...,kn) = > 0. 

Denote now by di, i = 1, . . . ,n the diagonal entries of and let us shown next that 
di > ^ for all n > 1, for alH = 1, . . . , n. To see this, for each i G {!,..., n}, denote by 
Ci the z-th unit vector of the basis of the Euclidean space and consider the coordinate 
change matrix T = [e* ^2 - ei_i ei e,+i - e„]^ go that Tj = = and so that TiBTi 
exchanges the 1-st and the i-th rows and columns of B. Then it is readily seen that TiBTi 
is positive definite (by positive definiteness of 5, as established in the proof of Lemma |3] 
in Section I5.2p and that the (1, 1) element of TiB~^Ti = {T.iBT.j)~^ corresponds to d^. 

\ b b"^ 1 

Partition the matrix as TjBTj = ^ ' , where 6,- is the i-the diagonal element of B 

satisfying by construction bi < 2. Since -64,22 > 0, the following holds by applying the 
matrix inversion formula to the (1,1) element of TiB~^Ti = {TiBTi)~^: 



di = [bi 



C2 5.iM"'>^r'>2-'- 



(10) 



As a next step, considering that ctj, i = 1, . . . ,n are distinct positive even integers, we 
have that 

n ^ n ^ ^ n 



^ a? - 1 - ^ (2j)2 _ 1 ~ 2 ^ 2? - 1 2? + 1 



1 



1 



^ - 2j - 1 2j + 
2 V 2n + 1 / - 2 



;ii) 



Finally, keeping in mind that the trace of KB ^ is the sum of its eigenvalues and consid- 
ering equality ([7]) we can combine ( ITU]) and ffTT]) to get 



^ n n n 

oE^^<E^A = Tr(i^i?^^) = $^ 



j=i 3 



1 
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which, bearing in mind that kj > 0, for all j G {1, . . . , n}, implies, for each i G {1, . . . , n}, 
ki < ^"=1 kj < 1 as to be proven. □ 

Remark 3 (Sensitivity analysis) Sensitivity of the solution obtained by numerical tech- 
niques (or also by symbolic techniques, by truncating rational numbers to floating point 
numbers) can be assessed a posteriori. Indeed, the eigenvectors of a matrix encode the 
sensitivity of the eigenvalues to (unstructured) uncertainty affecting the entries of the ma- 
trix itself. In particular, sensitivity can be and will be assessed numerically in Section 13.41 
using two main methods: a first one providing local only information by using the condi- 
tion number of simple eigenvalues given in [TTf §7.2.2] and a second one providing an idea 
about the effect of large perturbations by graphically displaying the pseudo-spectrum of 
the matrix (namely the sublevel sets of the perturbed eigenvalues under norm bounded 
unstructured perturbations) according to [32]. As shown in Section [3^ these tools provide 
useful insight about the features of different solutions leading to interesting interpretations 
to some of the observations in [B]. o 



3 Symbolic solution 

In this section we use techniques from real algebraic geometry to solve symbolically the 
inverse eigenvalue problem ([2]). We focus on system (issued from (jlj)) whose polyno- 
mials Pi conveniently inherit some sparsity of the tridiagonal matrix B, in the sense that 
many monomials of the fi variables are zero. In contrast, the polynomials qi in system 
iQ (issued from ([7])) are less sparse since matrix is dense. 

3.1 Assumptions 

Solving problem ([6]) amounts to finding a real n-dimensional solution to a system of n 
polynomial equations with rational coefficients. Let us denote this system as 

Pi(/i, •••,/«) 

P2(/l, •••,/«) 
Pnifl, ■ ■ ■ , fn) 

where pi are given scalar-valued multivariate polynomials and fi are scalars to be found, 
i = 1, . . . ,n. 

Throughout the paper we make the following standing assumptions. 

Assumption 1 There exists a finite number of complex solutions to the system of poly- 
nomial equations [W\) . 

Note that Assumption [T] readily implies that there is a finite number of real solutions to 
the system of polynomial equations (|T2|) (equivalently, ([6])). 








(12) 
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Assumption 2 There exists at least one real solution to the system of polynomial equa- 
tions (Q^j- 

In the conclusions section we state open problems in connection with these assumptions. 
3.2 Real algebraic geometry 

For an elementary tutorial account of real algebraic geometry, please refer to [Tl] . In this 
paragraph we survey only a few essential ideas, with a focus on explicit algorithms. 

Finding a solution to polynomial system (|T2|) amounts to finding a point in the set 

r :={/gC" : pi(/) = 0,...,p„(/) = 0}. 

Set 1^ is a subset of C" called a (complex) algebraic variety because it is the vanishing locus 
of a finite number of polynomials. To the geometric object Y corresponds an algebraic 
object: 

^:={ai(/)Pi(/) + --- + a„(/K(/) : ai, . . . , a, G Q[/]} 

which is a subset of Q[/] called an algebraic ideal. Elements in are obtained by taking 
linear combinations (with polynomial coefficients) of polynomials Pi. We say that ideal 
is generated by pi, . . . ,p„, and we say that pi, . . . ,p„ is a presentation of ^ . Although 
has an infinite number of elements, it follows from a fundamental theorem of Hilbert 
[TTl Theorem 4, §2.5] that it has always a finite presentation, and this allows to be 
handled by a computer. Note that every polynomial in vanishes at points / G and 
we say that Y is the variety associated to the ideal ^ . 

A key idea of algebraic geometry consists in obtaining useful information on Y from a 
suitable presentation of ^ . By taking finitely many linear combinations of polynomials pj, 
we will generate another equivalent system of polynomials which generates the same ideal 
but with another presentation, and which is associated with the same variety 'Y . In 
particular, when is a discrete set, i.e. when Assumption [1] is satisfied, this presentation 
should allow to compute the solutions easily. A useful presentation is a Grobner basis 
[TTl Section 2]. To obtain such a basis, we can devise an algorithm using only linear 
algebra and performing a series of multivariate polynomial divisions. These divisions can 
be carried out provided one defines a suitable ordering on the set of monomials of variables 
fi, . . . , fn- For our purpose of computing the solutions, a useful ordering is the graded 
reverse lexicographic (grevlex) order, see [HI Definition 6 of Section 2.2]. Once a Grobner 
basis is available, we can compute a rational univariate representation (RUR) 



n(/o) 

-(/°) (13) 

rnjfo) 

r{fo) 

where r is a univariate polynomial in a variable /o, and ri, . . . , r„ are univariate polyno- 
mials of degree less than the degree of r. Variable fo is called the separating variable, and 
it is a linear combination of the original variables /i, . . . , /„. 



r{fo) = 

h = 

fn 
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Proposition 1 For system [T^) . a Grobner basis always exists. Moreover, all but a finite 
number linear combinations of the variables fi, i = l,...,n are separating variables. 
Finally, once the separating variable is chosen, the RUR exists and is unique. 

Proof: The existence of a Grobner basis follows from [TT| Corollary 6, §2.5]. The rest 
of the proposition can be proven using [TB], [271 Theorem 3.1] or [21 Proposition 12.16]. 
See also [3T| Proposition 2.3] and the discussion just after, which explains the connection 
between RUR and Grobner basis. □ 

Once a RUR is available, enumerating all the real solutions amounts to computing all the 
real roots of the univariate polynomial r in (fT3|) . and evaluating the rational functions 
Vi/r in ( IT3|) at those points. By an appropriate use of Descartes' rule of signs and Sturm 
sequences, see e.g. [21 Section 2.2], an algorithm can be designed that isolates all real 
roots of the univariate polynomial in rational intervals of arbitrarily small width. 



3.3 Maple implementation of the algorithm 

To solve the polynomial system of equations (IT2]) we use special computing packages which 
can be called directly from inside a computing sheet of the Maple computation software. 

The first step computes a Grobner basis in grevlex ordering. This is achieved by means 
of the fgb_gbasis procedure of the FGB package, see [TH [15]. This implementation of 
Grobner basis computation is considered to be one of the most efficient available. To 
speed up linear algebra computations, the algorithm works in an integer ring with prime 
characteristic, for a sufficient large prime number found iteratively. Then, for finding the 
real solutions, a RUR is computed from the Grobner basis, see [27]. Finally, for isolating 
the real roots of the univariate polynomial, the procedure rs_isolate_gb in the f gbrs 
package is used. It returns small rational intervals (as small as one wants) within which 
the roots are guaranteed to be found. See [28] for further references on the method used. 

The whole process lies in a few lines of Maple code. We first generate the polynomial 
system (corresponding to ([6|)), denoted by p in the following Maple sheet: 

with(LinearAlgebra) : with(PolynomialTools) : 
n:=4:B:=Matrix(n) : 
for i from 1 to n-1 do 

B(i,i):=2: B(i , i+1) : =-1 : B(i+1 , i) : =-1 : 
end do : 

B(n,n) :=(n+l)/n; 

F :=Matrix(n, Vector (n,syinbol=f) , shape=diagonal) : 
d:=product(x-((2*j)"2-l) , j=l. .n) : 
p : =CoefficientList (collect ( 
charpoly(B.F,x)-d,x) ,x) ; 

We then use the FGb package to transform the system p in a new algebraic system refer- 
enced with the name GB: 
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with(FGb) : 

fv: = [seq(f [i] ,i=l. .n)] : 
GB : =f gb_gbasis (p , , f v , [] ) : 

Finally, the solutions are computed as follows: 

with(fgbrs) : 
rs_isolate_gb(GB,fv) ; 

As an example, in the case n = 4 the univariate polynomial r(/o) that enters the RUR 
(|T3|) has degree 24, and it parametrizes all the 4! = 24 (complex) solutions. Out of 24, 
only 4 solutions are real. This polynomial, or rather its scaled version 

fo ^ sign(r(/o)) |r(/o)|^, 

which allows to more easily visualize its zero crossings, is represented in Figure |2l It is 
computed as follows 

with(Groebner) : 

R: =RationalUnivariateRepresentation(GB ,x) : 
r:=R[l] ; 
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Figure 2: Case n = 4: rescaled graph of the univariate polynomial whose 4 real roots 
parametrize the 4 real solutions of the polynomial equations (E]) . 
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In practice, the computation of a Grobner basis becomes very hard as n, the number of 
variables and equations, increases. For the specific system of equations (E]) (or, equiva- 
lently, f lT2]) ) associated to the formulation (jl]) of the Marx generator design, we observe 
that for 2 < n < 6 the univariate polynomial r{fo) of the RUR has degree n\. Moreover, 
the size of the (integer) coefficients of this polynomial are very large. If one compares 
the degree of the polynomial, 720 in the case n = 6, and the number of real roots, 12 in 
this case, it is clear that the computation complexity is due to the very large number of 
complex roots, which are of no interest for our engineering problem. 

3.4 Numerical results and sensitivity analysis 

The simplest possible selection of parameters in the assignment problem ([2]) is to select 
them as the smallest possible set of distinct even and positive integers, namely a, = 2i for 
alH G {1, . . . , n}. This selection gives rise to the smallest possible coefficient list in ([6]) and 
leads to a set of solutions to Problem [1] having the smallest possible maximal frequency of 
the natural resonant modes of the circuit. This was also the preferred solution addressed 
in [5^ E]. In particular, for this specific selection of the resonant frequencies, a set of 10 
solutions for the case n = 6 has been given in [6l Table II]. The advantage of the formal 
approach of this section allows to find the complete set of solutions, amounting to 12 (2 
solutions were missed in [6l Table II] which used bounded random sampling followed by 
numerical optimization). Another advantage of our results is that the sets of solutions 
reported here were computed in a few minutes using the Maple code reported in the 
previous section. Table [1] reports the numerical values computed by the analytic algorithm 
proposed in this section. Note that the displayed values correspond to v?Ci/c = v^f^^. 
We choose to represent these values to allow for an easier comparison with the results in 
Table II] and because they are better numerically conditioned. Note that solutions 
number 1 and 2 for the case n = 6 were not reported in [51 Table II]. 

According to the observations reported in Remark [3] we can characterize the sensitivity 
of each solution obtained from the proposed symbolic solution method by looking at the 
condition of each eigenvalue of matrix BF in ([2]). This method can be applied because all 
eigenvalues of BF are distinct (therefore simple) by assumption. In particular, for each 
eigenvalue of BF its condition number corresponds to where w and v have unit 

norm and are respectively the left and right eigenvectors associated with that eigenvalue. 
The results of the sensitivity analysis is represented by the maximum condition among 
all eigenvalues of matrix BF and is shown in the last column of Table [T] for each one of 
the computed solutions. 

Inspecting the different sensitivities it appears that the last solution for each one of the 
analyzed cases corresponds to the least sensitive one, namely the one that is expected 
to be more robust. Interestingly, this solution corresponds to the solution qualitatively 
characterized in [6] as the "regular" solution. Indeed, when looking at the time responses 
of the Marx generator network designed with these parameters, one experiences little 
dependence on higher order harmonics and some suitable monotone evolution of certain 
voltages in the circuit (see [6l Fig. 10] for an example of this). Another peculiar feature of 
the "regular" solutions corresponding to the last solution for each n in Table [1] is that the 
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1 
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n = 2 














1 


1.50213 


0.47340 










1.11020 


2 


0.63120 


1.12660 










1.02661 


n = 3 














1 


1.49303 


1.49229 


0.41548 








1.15573 


2 


0.84408 


0.77662 


1.41217 








1.03869 


n = 4 














1 


1.71070 


1.29555 


1.54667 


0.38529 






1.18493 


2 


1.62637 


0.62519 


1.73498 


0.74862 






1.09175 


3 


1.06181 


2.10211 


0.66491 


0.89099 






1.12097 


4 


1.13210 


0.78731 


0.92450 


1.60306 






1.04400 
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cond 
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C 




n 


= 5 














1 


2.04567 


1.23900 


1.35694 


1.57111 


0.36796 




1.20185 


2 


2.14782 


0.63778 


1.24610 


1.70028 


0.68506 




1.10925 


3 


0.99720 


1.69936 


1.57266 


0.64267 


1.16088 




1.05580 


4 


1.47480 


0.86342 


0.84481 


1.07344 


1.72179 




1.04479 


n 


= 6 














1 


2.49095 


1.25588 


1.20240 


1.49359 


1.53290 


0.35987 


1.20655 


2 


1.92537 


1.79971 


1.80083 


0.90696 


1.45858 


0.37545 


1.19542 


3 


1.67555 


1.98118 


2.05786 


0.66667 


1.30433 


0.52174 


1.11204 


4 


2.65073 


1.01602 


0.68610 


1.82261 


1.42150 


0.64738 


1.15057 


5 


1.34706 


2.18478 


0.92044 


1.84208 


0.65432 


0.94921 


1.09706 


6 


1.95229 


0.93587 


1.53272 


0.63062 


1.76898 


0.99206 


1.08439 


7 


2.46541 


0.73028 


0.99423 


0.93809 


1.85839 


0.99314 


1.12232 


8 


1.79820 


0.94167 


1.74742 


0.62528 


1.59040 


1.05327 


1.08837 


9 


1.43355 


1.89698 


0.60976 


1.73302 


1.02388 


1.05334 


1.10087 


10 


1.50458 


1.00778 


2.04896 


1.05339 


0.68756 


1.37732 


1.06859 


11 


1.38734 


1.13092 


1.48392 


1.32482 


0.63764 


1.57578 


1.06168 


12 


1.87892 


0.96056 


0.85587 


0.91518 


1.23619 


1.77345 


1.05917 



Table 1: The solutions to the Marx design problem computed using Grobner basis methods 
and their conditioning. 

interpolated values of n^Cj/c form a convex function of i, namely (since n^/c is constant) 
one has cj < ^^+'^~^^^''^ for all j = 2, ... ,n — 1 (see also the red curve in [6l Fig. 8] 
corresponding to the last solution for = 6 in Tabled]). Moreover, at least up to = 6, 
numerical evidence reveals that there only exists one such solution. Due to its desirable 
features both in terms of numerical robustness and of desirable time evolution of the 
arising circuit (as reported in |6]), we will be imposing this constraint on the numerical 
optimization described in the next section, to be able to isolate that specific solution for 
the case n > 6 (or all of such specific solutions, if more than one of them exists). 

A final comment should be made about the use of the pseudospectra of matrix BF as a 
graphical tool to assess the sensitivity of each solution of Tabled] to larger perturbations. 
According to the results in [32] and the corresponding Matlab tool Pseudospectra GUI 
available in the EigTool package [32], given e > 0, the associated pseudospectrum of a 
matrix A corresponds to the following region in the complex plane: 

{s G C: det(sJ„ - A - A) =0, for some A G C'^'''', || A|| < e} 

(see Figure [3] for some examples). Surprisingly, the graphical aspect of the pseudospectra 



13 




Figure 3: The pseudospectra of the state transition matrix A for the four solutions cor- 
responding to 77, = 5. 

appear indistinguishable among all solutions for fixed n. Nevertheless, interesting results 
are obtained when applying this analysis to the state transition matrix A associated to the 
state-space model of the circuit (for fixed values of i and c) whose expression is given in 
f|T8|) . in Section Ism An example of the type of pseudospectra obtained for the case n = 5 
is represented in Figure [3] corresponding to the case c = £ = 1 . The figure shows that 
once again the "regular" solution is associated with the least sensitive scenario. Visually, 
this corresponds to the tightest pseudospectra (see the rightmost case in Figure E]) which 
can be best appreciated by inspecting the largest level set of the figure, corresponding to 
the selection e = lO*^'^ according to the legend to the right of the figure. The other values 
of n lead to similar results. 



4 Numerical solution 

In this section we use convex optimization techniques to find numerically the real solutions 
of our polynomial system of equation. As compared to the previous section, we focus here 
on the alternative formulation because for this formulation, according to Lemma [21 all 
the real solutions satisfy \ki\ < 1, i = 1,2, ... ,n. As explained e.g. in [20], for numerical 
reasons it is very important that the problem unknowns are scaled down to the unit 
interval. 



4.1 Problem formulation 



A numerical approach solution to the inverse eigenvalue problem presented in Section 12.31 
consists in formulating it first as a nonconvex polynomial optimization problem: 

q* = miufc qo{k) 

s.t. keje^ ^ ' 

where the objective function go G Q[/c] is a given polynomial of the vector of indeterminates 
k G M"', and, based on ([9]) and on the discussion about "regular solutions" in Section [331 
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the feasibility set 

= {keW : qi{k) = 0, z = 1, . . . , n, gj{k) := kj - 2kj+i + kj+2 > 0, j = 1, . . . , n - 2} 

is the real algebraic variety corresponding to the zero locus of the ideal studied in 
Section [3] intersected with the polyhedron modeling the convexity constraints introduced 
in Section [331 

A typical objective function in problem 0141) can be the positive definite convex quadratic 
form 

n 

q^{k) = Y,ih-k,f = k^ 

SO that capacitors Cj = ki/c are as identical as possible, but we can also consider other 
relevant objective functions, not necessarily quadratic, definite in sign or convex. 

Optimization problem fll4l) is finite-dimensional, algebraic, but nonconvex since the fea- 
sibility set ^ is disconnected, as a union of isolated points (because of Assumption [1]). 
Local optimization techniques based on nonlinear programming are likely to face troubles 
with such sets. 

Function go is continuous and we optimize it over which is compact, since by Assump- 
tions [U and [2] there is at least one real solution and at most a finite number of isolated 
real solutions. It follows that optimization problem f lT^ has at least one solution. 

Since go is not necessarily convex and is disconnected, we do not expect optimization 
problem (H^ to have a unique global minimizer. However, we expect the number of global 
minimizers to be significant smaller than the cardinality of set J^. 

4.2 Optimization method using Gloptipoly 

Our optimization method is based on an idea first described in [22] which consists in 
reformulating a nonconvex global optimization problem with polynomial data (i.e. min- 
imization of a polynomial objective function subject to polynomial inequalities and/or 
equations) as an equivalent convex linear programming (LP) problem over probability 
measures. Instead of optimizing over a vector in a finite-dimensional Euclidean space, we 
optimize over a probability measure in an infinite-dimensional Banach space. The mea- 
sure is supported on the feasibility set of the optimization problem, which is algebraic in 
our case, and we require in addition that the set is bounded (which is true in our case, by 
assumption). More concretely, a probability measure is understood as a linear functional 
acting on the space of continuous functions, and we manipulate a measure through its 
moments, which are images of monomials (which are dense w.r.t. the supremum norm 
in the space of continuous functions with compact support). Using results on functional 
analysis and real algebraic geometry, and under some mild assumption on the compact 
support, a sequence of real numbers are moments of a probability measure if they be- 
long to an appropriate affine section of the cone of positive semidefinite linear operators, 
an infinite-dimensional convex set. We then construct a hierarchy of finite-dimensional 
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truncations of this convex set, namely affine sections of the cone of positive semi defi- 
nite matrices of fixed size. As a result, solving an LP in the set of probability measures 
with compact semi- algebraic support boils down to solving a hierarchy of semidefinite 
programming (SDP) problems, also called linear matrix inequalities (LMIs). 

When there is a finite number of global optimizers, the approach is guaranteed to converge 
in a finite number of steps, and the global optimizer(s) can be extracted with the help of 
numerical linear algebra, see [20]. We have then a numerical certificate of global optimality 
of the solution(s). This approach has been successfully applied to solve globally various 
polynomial optimization problems, see [23] and [23] for general overviews of results and 
applications. For applications in systems control, the reader is referred to the survey [T^ . 
For finding real solutions of systems of polynomial equations and real radical ideals, the 
approach has been comprehensively studied in [23]. 

More explicitly, we now describe our approach to the numerical solution of problem (1141) . 
We consider compact set C and we denote by ^(J^) the Banach space of Borel 
measures supported on J^. These are nonnegative functions from the Borel sigma-algebra 
of subsets of to the real line M. Given a measure /x G ^( J^) we define its moment of 
order a G N" as the real number 

ya= x'^fiidx) e R (15) 

where we use the multi-index notation for monomials, i.e. = YYi=ix'^'. We define the 
infinite-dimensional vector y = {ya}aeN" as the sequence of moments of /i. Note that 
yo = J fJ' = fJ'i'^) = 1 whenever fi G ^{J^) is a probability measure. Moreover, if for 
each k G J^, \ki\ < 1 for all i (this is what we establish in Lemma [2]), then \ya\ < 1 for 
all a G N™. Conversely, for larger sets J^, the variable ya may grow very large and this 
is not convenient for numerical reasons. This aspect has been pointed out in [20] and 
motivates Lemma |21 

Given a sequence y, we define the Riesz linear functional iy : — M which acts 

on polynomials 7i{x) = J2a 

TTfyX clS follows: £y{7i{x)) = J^a'^aVa- If scqucnce y has a 
representing measure /i, integration of polynomial 7i{x) w.r.t. /i is obtained by applying 
the Riesz functional iy on 7r(x), since iy{7i{x)) = J 7r{x)fi{dx) = j Ylia^^iX" ^{dx) = 

If we apply the Riesz functional on the square of a polynomial 7r(x) of degree d, then we 
obtain a form which is quadratic in the coefficient vector tt = {'n:Q}\a\<d and which we 
denote 

iy{T:\x))=T:^Md{y)Ti 

where M^iy) is a symmetric matrix which is linear in y, called the moment matrix of 
order d. Rows and columns in this matrix are indexed by vectors a G N" and /3 G N", 
and inspection reveals that indeed the entry (a,/3) in matrix Md{y) is the moment ya+i3- 
Given a polynomial x{^) we let 

4(7r2(x)x(x)) = 7r^Mrf(x,y)vr 

^The notation ai is used in (jlSp and the remaining derivations in this section, for consistency with 
the notation used in |22) and references therein. However, they should not be confused with the scalars 
ai used in Theorem [TJ 
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where M^lx, y) is a symmetric matrix which is hnear in y and hnear in coefficients of 
x{x), called the localizing matrix of order d w.r.t. The localizing matrix is a linear 

combination of moment matrices, in the sense that entry (a,/?) in Md{x,y) is equal to 

4.3 Application to the eigenvalue assignment problem 

Based on the optimization method presented in the previous section, we now formulate the 
polynomial optimization problem f lT^ as a hierarchy of finite-dimensional LMI problems 
with the help of the Matlab interface GloptiPoly 3 [21]. Then, we use public-domain im- 
plementations of primal-dual interior point algorithms to solve the LMI problems. These 
algorithms rely on a suitable logarithmic barrier function for the SDP cone, and they 
proceed by iteratively reducing the duality gap between the primal problem and its dual, 
which is also an LMI problem. Each iteration consists of solving a Newton linear system 
of equations, involving the gradient and the Hessian of a Lagrangian built from the bar- 
rier function. Most of the computational burden comes from the construction and the 
storage of the Hessian matrix, and problem sparsity can be largely exploited at this stage. 
For more information on SDP and related optimization methods, see e.g. [1]. For our 
numerical examples we have been using the SDP solver SeDuMi 1.3 |30j . 

More specifically, let di = \ '^^^ ^' ] , i = 0, 1, . . . , n, and consider the optimization problem 

1*d = inf?; ^y{qo) 
s.t. yo = l 

Md{y) h (16) 

Md-diiQi,y) = i = l,...,n 

Md-i{gj,y)hO, j = l,...,n-2 

for d > max{(ij}j=o,i,....n, where >z stands for positive semidefinite. 

In the above problem, the unknown is the truncated sequence y of moments of degree up 
to 2d, and the constraints are convex linear matrix inequalities (LMI) in y. Problem f|T6|l 
is called the LMI relaxation of order d of problem f|T^ . It can be proved that the infimum 
in LMI problem ([TB|) is attained, as stated in the next proposition which is proven in [251 
Theorem 6.1]. 

Proposition 2 The optimal values of / f7^) satisfy < Qd+i- Moreover, under Assump- 
tionsUl andl^ there exists a finite d* such that = q* for all d> d* , where q* is the 
optimal value of [I4\ )- 

We should remark that while LMI solutions are not very accurate, they can be obtained 
cheaply (at least for these examples) since we do not need to enumerate all real solutions, 
only the optimal one(s). Moreover, if a computed solution is not deemed accurate enough, 
it can be refined locally afterwards by Newton's method if required. 

Remark 4 Proposition [2] states that solving nonconvex polynomial optimization f|T4l) is 
equivalent to solving convex LMI problem f|T6l) for a sufficiently large relaxation order d* 
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if the feasibility set has zero dimension and is nonempty . Unfortunately, it is not possible 
to give a priori useful (lower or upper) bounds on d*, and the strategy followed in [20] is 
to detect global optimality of an LMI relaxation by inspecting the rank of the moment 
matrix and then extract the globally optimal solutions by linear algebra, see also [21] and 
[23| Algorithm 6.1]. For all the Marx generator design problems that we considered, we 
observed that the global optimum is certified at the smallest possible LMI relaxation, i.e., 
in Proposition 121 

d* 





degg. 






max 






i=0,l,...,n 


2 




2 



4.4 Maple/Matlab code example for n = 3 

First we generate the polynomials qi{k), i = 1,. . . ,n corresponding to Q and defining 
the feasibility set in problem (|T4l) . We use the following Maple code, where we select 
ai = 2i, i = 1, . . . ,n just like in Section 13.41 

with(Linear Algebra) : with(PolynomialTools) : 
ii:=3:B:=Matrix(n) :for i from 1 to n-1 do 

B(i,i):=2: B(i , i+1) : =-1 : B(i+1 , i) : =-1 : 
end do: B(n,n) :=(n+l)/n; 
K : =Matrix (n , Vector (n , symbol=k;) , 

shape=diagonal) : 
p:=product(x-l/((2*j)~2-l) , j=l. .n) : 
q: =CoefficientList (collect ( 

charpoly(MatrixInverse(B) .F,x)-p,x) ,x) ; 

For n = 3 this code generates the following polynomials 

qi{k) = -f + |/ci + |/C2 + 1^3 
q2{k) = -^ + ^kik2 + kik3 + k2k3 
qsik) = -j^ + ^kik2k3. 

These polynomials are then converted into Matlab format, and we use the following Glop- 
tiPoly code for inputing problem f|T^ and solving the smallest possible LMI relaxation, 
i.e. d = [|] = 2 in problem f lTB]) (note also the inequality constraints in P corresponding 
to the convexity requirement for the "regular" solution): 

mpol k 3 

K = [-3/7+5/6*k(l)+4/3*k(2)+3/2*k(3) 

-53/1575+2/3*k(l)*k(2)+k(l)*k(3)+k(2)*k(3) 

-l/1575+l/2*k(l)*k(2)*k(3)] ; 
obj =0; 

for i = l:length(k) 
for j = l:length(k) 

obj = obj+(k(i)-k(j))-2; 
end 
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end 

P = msdp(min(obj) ,K==0,k(l)-2*k(2)+k(3)>=0) ; 
[stat,obj] = msol(P); 
double (k) 

As pointed out in Remark HJ already with the smallest LMI relaxation (i = 2, we obtain 
a certificate of global optimality, and a unique global minimizer (truncated at 5 digits): 
{ku k2, ks) = (9.3786 ■ IQ-^, 8.6296 • IQ-^, 1.5690 ■ 10"^) that corresponds to the second 
solution in the second block of Table [D (the "regular" one as expected). 

4.5 Numerical results and sensitivity analysis 

We applied this numerical approach on our examples, with different objective functions, 
and the overall conclusion is that the cases n = 2, 3, 4, 5 are solved very easily (in a few 
seconds) but the solution (obtained with SeDuMi) is not very accurate. The case n = 6 
is solved in a few minutes, and the case n = 7 is significantly harder: it takes a few hours 
to be solved and it provides the following "regular" solution (in the sense introduced in 
Section [331): 



n'' — 

c 


2C2 

n — 

c 


2C3 

c 


2C4 

n — 

c 


2C5 

n — 

c 


2C6 

c 


2C7 

n — 

c 


2.07061 


1.05669 


1.04940 


1.05715 


1.06861 


1.08449 


1.85298 



For this solution, we can compute the sensitivity level using the same algorithm used in 
the last column of Table [1] and we obtain 1.0502. 

Finally, solving the case n = 8 takes approximately 15 hours and leads to the following 
set of parameters. The sensitivity of this solution corresponds to 1.0617. 



2C1 
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c 


2^3 
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2C4 
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2C5 

n — 

c 


2^6 

n — 

c 


2C7 
C 


2^8 

n — 

C 


2.39407 


1.17326 


1.12475 


1.11221 


1.10440 


1.09960 


1.09985 


1.87282 



5 Circuit description and proofs 

We carry out the proof of Theorem [T] by first providing a mathematical description of 
the circuit (Section 15. ip . then proving a useful resonance result (Section 15.21) and then 
proving the theorem (Section 15. 3p . 

5.1 Circuit description 

Following similar approaches to those adopted in [1], it is possible to derive a state-space 
representation of the Marx generator of Figure [1] using 3n + 2 state variables comprising 
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Figure 4: Sign conventions in the models of the Marx generator circuit. 



2n + 1 voltages across the 2n + 1 circuit capacitors, and n + 1 currents flowing in the n + 1 
inductors. 

In particular, using the sign convention and the component values depicted in Figure S] 
(see also Figured]), we can define the state variable as 

T 



and the linear dynamics of the circuit correspond to the equations 

cvk = ik, k = 1,. . n, cvn+i = 
Lii = I'd - t'l, nLiL = Vcn - nvn+i, 
Lik = Vck - Vc{k-i) -Vk, /c = 2, . . . , n, 
4+1 



(17) 



CkVck 



^5 k — 1, . . . , Ti 



which can be written in compact form using the following linear state-space model 



X = Ax = 


0(2n+l)x(2n+l) 


1 T 

-cj:f"^j_i 


X, 


(1 




It lyT 7-1 


0(n+l)x(n+l) 







where Oqxq G W^^'^ is a square matrix of zeros, Iq is the identity matrix of size q 

Jn+i = diag(l,...,l,n) G m("+i)x("+i), = diag(l, . . . , 1, -1) G 

ri -1 ■■■ On 
1 -1 ■■• 



agonal matrices, S 
Theorem [H 



LO 



1 -IJ 



are both di- 

G M"x("+^) and F is defined in the statement of 



An alternative route can be taken by describing the loop equations of the circuit in the 
transformed Laplace domain. In particular, using the fact that all the initial currents 
are zero and that the initial voltages across the storage capacitors is Vo while the initial 
voltages across the parasitic capacitors is zero, we get the following Kirchhoff's voltage 
law equations for the n + 1 loops of the circuit: 



^Ji + s£h + ^ (/i - h) = ^ 

5C SCl ^ ^' S 

(Ik - Ik-l) + j-Jk + silk 



'fc+lj 



2,...n, 



(19) 



n+l 



Snlln+l + 77 /n+l = 0, 



where /^(s), = 1, . . . n are the Laplace transforms of the loop currents («i(t), . . . , i„(t), — «„+ 
(the last loop current has been changed in sign for convenience of notation) and the de- 
pendence on the Laplace variable s is omitted for compactness of the equations. 
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The above equations can be written in matrix form as 



G(s)I(s) = h-, 
s 

where G{s) e M("+i)x("+i), /(s) g M"+i and b G are defined as I{s) := 
and 



G{s) 



di 







1 

SCl 
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sc„ 



" 




"l" 








1 





b := 




]_ 




1 


dn+l _ 








^ + si + ^, 

SC SCl 

i + + + 2<k<n 

SC ■5Cfc_l SCk ' 



k = 1 
2 < A; 
/c = n + 1. 



(20) 



(21a) 



(21b) 



Consider now the fact that 1^+1(0) = 0. Then in the Laplace domain we have Vl{s) = 
— — /„_i_i(s) (recall that J„+i is the Laplace transform of —in+i), which can be combined 
with to get 



77 

K+i(s) = -^EG{s)-'bVo =: W{s)vo, 



(22) 



where the output matrix = [O ■ ■ ■ l] ensures EI{s) = In+i{s). The scalar transfer 
function W{s) is the transfer function from the specific initial condition characterized in 
Problem [T] to the output corresponding to the voltage across the load capacitor. 



5.2 A sufficient resonance condition 

In this section we establish a preliminary result that will be useful for the proof of The- 
orem [1] and which are based on the description fl22|) of the relation between the initial 
voltage fo and the voltage Vn+i{s) across the last capacitor. 

Lemma 3 Consider the Marx circuit in Figure [I] and the matrices B in ^ and F = 
diag(c/ci, . . . , c/cn). The matrix I + BF has n real positive eigenvalues. Moreover, denot- 
ing by ai, . . . ,an such n real positive eigenvalues and fixing oq = 1, the transfer function 
W{s) in has at most 2n + 3 poles at s = and s = ±j-y/|^, k = 0, . . . ,n, before 
possible cancellations. 

Proof: Consider the matrix G{s) in f l2T|) and note that it can be written as G{s) = 
j^Gp{s), with Gp{s) being a polynomial matrix defined as: 

Gp{s) = (1 + s^ic)Jn+i + S^FE, (23) 
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with Jn+i, S defined below f|T8|) and F is defined in the statement of the lemma. To study 
G{s)~^ in (1221) . let us first compute the determinant of Ap = det(Gj,(s)), using the form 
fl2^ and using Sylvester's determinant theorem stating that det(/— X-'^F) = det(/— FX"^) 
for any two matrices X, Y having the same dimensions: 



Ap = det((l + s^ic)Jn+i + S^FS) = det (4+1 + S^F ) det((l + s^ic)Jn+i) (24) 

n(l + s^icr^' det (^In + = ^(1 + s'ic) det ((1 + s'ic)!^ + S J-_;,S^F) 

n(l + s'^ic) det{s'^icln - {-In - BF)), 



where Ik is the k-th dimensional identity matrix, and it can be easily verified that B = 
EJ~^':^E^ corresponds to the expression in ([1]). 

First we establish that I+BF has n positive eigenvalues. This follows from the coordinate 
change T = ^/F transforming the matrix into I + \fF B\/F, which is a symmetric positive 
definite matrix because B = TiJ^^^H'^ is positive semi-definite by positive definiteness of 
Jn+i- More specifically, it is easy to see that B is actually positive definite by decomposing 
S = [s -e„+i], where S is an upper triangular matrix and e„+i = [0 ■■■ 1]"^, then B = 
SS"^ + -^e„+ie^_,_^, where the first element is positive definite and the second one is positive 
semi-definite. 

Consider now the poles of W{s) in f l22]) and note that W{s) = —-EGp{s)~^b, where 
Gp{s)^^ = A~^Adj(Gp(s)), and Adj(Gp(s)) is the adjugate matrix of Gp{s) which is 
polynomial because Gp{s) is polynomial itself from (l23i) . Therefore the poles of W{s) are 
at most (before possible cancellations) the pole at the origin and the roots of Ap which, 
from (1241) appears to be a biquadratic polynomial of degree 2n -|- 2, thus proving that 
W{s) has at most 2n + 3 poles. 



As for the position of these poles, note that the first term in the last equality of 
becomes zero for s'^ic = —1 = —a^, while the second term becomes zero when s'^ic 
coincides with the eigenvalues of — (/ + BF), namely s'^ic = — 1 = —ak, k = 1, . . . ,n, 
which implies that the roots of Ap are given by s = ±jy^, k = 0, . . . ,n, plus the pole 
at the origin, as to be proven. □ 



5.3 Proof of Theorem [T] 

Consider the circuit of Figure [1] and its state space equations corresponding to f|T8|) with 
the state x given in (fT7|) . The state-space equations involve 3n + 2 eigenvalues (and 
eigenvectors), corresponding to the 2r2 + 2 resonant modes characterized in Lemma [3] and 
n eigenvalues at zero, which can be easily determined by setting all currents to zero and 
imposing a constant (integer, for convenience) voltage at each one of the Vck nodes (refer 
to Figure [1] as well) . 
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In particular, the following n eigenvectors correspond to the n eigenvalues at zero: 
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1 

Nn 

On+1 



(25) 



where Oq G M'^ is a vector of zeros, Iq G M'' is a vector whose entries are all 1, and Ng G M'' 
is a vector whose entries are increasing natural numbers from 1 to q. Indeed it is easy to 
check that AwQk = 0, for all /c = 1, . . . , n, where A is defined in fll8l) . 

While the n eigenvectors in f l25p are a basis for the n-dimensional eigenspace corresponding 
to all the eigenvalues of A, the main resonant mo d^ of the circuit, corresponding to the 

pair of eigenvalues ±j -y/^ is associated with an eigenspace whose basis can be written 
as: 
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indeed it can be easily verified that, using A as in 
cAw2 = Wi. 



we have iAwi = —W2 and 



Now, consider the initial condition characterized in Problem [T] and note that it can be 
written as: 



x{0) 
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On+1 
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On 



-Nn 
On+1 



^ {Wi + Won + W3) 



where Wi, Won are given in fl26l) and fl25l) . respectively, and W3 does not contain any 
component in the eigenspace of the main resonant mode (indeed its scalar product with 
both Wi and W2 is zero). 

Select now T = as in the statement of Theorem [H which corresponds to half of the 
period of the main resonant mode and note that, by linearity, the circuit response can be 
written as the sum of the three responses from the three initial conditions wi, won and 
W3, namely: 

" Or. 



x{T) 



Vo 



{-Wi + Won + W3) 



'n 
Vo 
On 
On+1 



(27) 



where the first terms exhibits a sign change from the initial condition because the response 
is a pure cosine and we sample it at half of its period; the second term is unchanged 



^Scc also the statement of Lemma |3] which proves, additionaUy, that this main resonant mode is always 
present in the cireuit, regardless of the selection of the parasitic capacitors c^, fc = 1, . . . , n. 
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because it lies in the eigenspace with zero eigenvalue; and the remaining term is also 
unchanged because its projection on the eigenspace of the main resonant mode is zero 
and by Lemma [3] (see also the reasoning at the beginning of this proof) all the remaining 
eigenvalues are either at zero (so that the corresponding response remains constant) or are 
pure imaginary with frequency being a multiple of This last property holds because 
from the statement of Theorem [H the integers a^, k = 1, . . . ,n are all even. And then, 
comparing ([2]) with the statement of Lemma [31 it appears that ak = ak, k = 1, ... ,n so 
that Ofc even for all /c > 1 implies (by Lemma [3]) that the frequency of all the remaining 
resonant modes is an even multiple of the main resonant mode. 

The proof is completed by noting that the state x{T) in (12 7p corresponds to having 
Vn+i{T) = Vo and all other voltages and currents at zero, which in turn implies vl{T) = 
nVn+i{T) = nvo, as to be proven. 

6 Conclusion and perspectives 

We proved that the design of an n-stage Marx generator electrical network exhibiting 
a desirable energy transfer boils down to a structured pole assignment. This can be in 
turn formulated as a structured system of n polynomial equations in n unknowns with 
rational coefficients. We have then illustrated a symbolic and a numerical approach to 
the computation of its solutions. By extrapolating from the analyzed cases, we conjecture 
that there is a finite number of complex, hence real solutions to this polynomial system. 
We also conjecture that there is at least one real solution. This motivates our Assump- 
tions [U and [21 The degrees of the polynomials are equal to 1, 2, . . . , n so that the Bezout 
bound on the number of complex solutions, as well as the mixed volume of the support 
polytopes of the polynomials [311 Section 3], both equal n\. The number of computed real 
solutions is however much less than this upper bound. It would be insightful to study 
the applicability of existing upper and lower bounds on the number of real solutions of 
systems of polynomial equations, see [22] for a recent survey. 

We solve the polynomial system of equations first with a symbolical method, namely 
Grobner bases and rational univariate representations. All real solutions are then obtained 
from the real roots of a univariate polynomial of degree n\. This univariate polynomial 
can be computed exactly (i.e. with rational coefficients) and its real roots are isolated in 
rational intervals at any given relative accuracy. Using state-of-the-art implementation 
of Grobner basis algorithms, we could solve the equations up to n = 6 stages routinely 
on a standard computer. Then we solved the same system of polynomial equations with 
a numerical method, using Lasserre's hierarchy of convex LMI relaxations for polynomial 
optimization. The advantage of this approach is that it is not necessary to represent or 
enumerate all n! complex solutions, and a particular real solution optimal with respect 
to a given polynomial objective function can be found quite easily up to n = 8. The 
accuracy of the computed solutions may not be very good, but this can be refined locally 
afterwards using, e.g. Newton's method. Future work may involve the use of alternative 
solution methods such as, e.g., homotopy or continuation algorithms [33l IT8| [3]. 
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